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Abstract 

We explore in some detail the hypothesis that the generation of a primordial lepton- 
antilepton asymmetry (Leptogenesis) early on in the history of the Universe is the 
root cause for the origin of matter. After explaining the theoretical conditions for 
producing a matter-antimatter asymmetry in the Universe we detail how, through 
sphaleron processes, it is possible to transmute a lepton asymmetry - or, more 
precisely, a (B-L)-asymmetry - into a baryon asymmetry. Because Leptogenesis 
depends in detail on properties of the neutrino spectrum, we review briefly existing 
experimental information on neutrinos as well as the seesaw mechanism, which 
offers a theoretical understanding of why neutrinos are so light. The bulk of the 
review is devoted to a discussion of thermal Leptogenesis and we show that for the 
neutrino spectrum suggested by oscillation experiments one obtains the observed 
value for the baryon to photon density ratio in the Universe, independently of any 
initial boundary conditions. In the latter part of the review we consider how well 
Leptogenesis fits with particle physics models of dark matter. Although axionic 
dark matter and Leptogenesis can be very naturally linked, there is a potential 
clash between Leptogenesis and models of super symmetric dark matter because the 
high temperature needed for Leptogenesis leads to an overproduction of gravitinos, 
which alter the standard predictions of Big Bang Nucleosynthesis. This problem 
can be resolved, but it constrains the supersymmetric spectrum at low energies and 
the nature of the lightest supersymmetric particle (LSP). Finally, as an illustration 
of possible other options for the origin of matter, we discuss the possibility that 
Leptogenesis may occur as a result of non-thermal processes. 
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1 Introduction 

Our understanding of the Universe has deepened considerably in the last 25 years, so 
much so that a standard cosmological model has emerged [1]. In this model, after the 
Big Bang, a period of inflationary expansion [2] ensued that effectively set the Universe's 
curvature to zero. After inflation the Universe's expansion continued, not in an exponen- 
tial fashion but with the rate of expansion being determined by which component of the 
Universe's energy density dominated the total energy density. 

In the present epoch this energy density is dominated by a, so-called, dark energy 
component whose negative pressure causes the Universe's expansion to accelerate [3]. 
Dark energy now accounts for approximately 70% of the total energy density, with the 
other 30% of the remaining energy density of the Universe's being dominated by some 
kind of non-luminous (dark) matter. In detail, the angular distribution of the temper- 
ature fluctuations of the microwave background radiation measured by the Wilkinson 
Microwave Anisotropy Probe (WMAP) collaboration [4] determines the various compo- 
nents of the ratio of the Universe's energy density now p to the critical energy density 
p c , Q = p j ' p c } The results are: f2dark energy = 0.73 ± 0.04; f2 ma tter = 0.27 ± 0.04; and 
Qb — 0.044±0.004. Here is the contribution of baryonic matter to f2 matter , confirming 
that about 85 % of fi matter is indeed contributed by dark matter. The contribution of 
neutrinos and photons is at a few per mil, or below, and is negligible. 

Although the standard cosmological model sketched above provides an accurate de- 
scription of the present Universe and its evolution, deep questions remain to be answered. 
What exactly constitutes the dark energy? Is it just a cosmological constant? But if that 
is so, why is the energy scale associated with the corresponding vacuum energy density 
[p cc = Eq] E ~ 2 x 10~ 3 eV] so small? Equally mysterious is the nature of the dark 
matter, although in this case there are at least some particle physics candidates that 
may be the source for this component of the Universe's energy density. 

A further mystery is associated with the observed baryon energy density. This num- 

x The critical density p c is the density that corresponds to a closed Universe now, p c — 3H% /8itGn ■ 
Here H a is the value of the Hubble parameter now. Inflation predicts that p a — p c , so that fi = 1. 
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ber can be used to infer the ratio of the number density of baryons to photons in the 
Universe, a quantity that is measured independently from the primordial nucleosynthesis 
of light elements. The WMAP results [4] are in agreement with the most recent nucle- 
osynthesis analysis of the primordial Deuterium abundance, but there are discrepancies 
with both the inferred 4 He and 7 Li values [5]. These latter values, however, may have an 
underestimated error [6]. Averaging the WMAP result only with that coming from the 
primordial abundance of Deuterium gives: 

-=?/b = 6.1±0.3x 1(T 10 . (1) 

n 7 

Why does this ratio have this value? 

In this review, we will principally try to address this last question which, as we shall 
see, is intimately related to the existence of a primordial matter-antimatter asymmetry. 
Nevertheless, we shall try, when germane, to connect our discussion with the broader 
issues of what constitutes dark energy and dark matter. 

There is good evidence that the Universe is mostly made up of matter, although it is 
possible that small amounts of antimatter exist [7]. However, antimatter certainly does 
not constitute one of the dominant components of the Universe's energy density. Indeed, 
as Cohen, de Rujula, and Glashow [8] have compellingly argued, if there were to exist 
large areas of antimatter in the Universe they could only be at a cosmic distance scale 
from us. Thus, along with the question of why %/n 7 has the value given in Eq. (1), 
there is a parallel question of why the Universe is predominantly composed of baryons 
rather than antibaryons. 

In fact, these two questions are interrelated. If the Universe had been matter- 
antimatter symmetric at temperatures of 0(1 GeV), as the Universe cools further and 
the inverse process 27 — > B + B becomes ineffective because of the Boltzmann factor, the 
number density of baryons and antibaryons relative to photons would have been reduced 
dramatically as a result of the annihilation process B + B — > 27. A straightforward 
calculation gives, in this case, [9]: 

n± = m „ 10 -is (2) 

Thus, in a symmetric Universe the question is really why observationally ns/n 7 is so 
large! 

It is very difficult to imagine processes at temperatures below a GeV that could 
enhance the ratio of the number density of baryons relative to that of photons much 
beyond the value this quantity attains when baryon-antibaryon annihilation occurs. 2 



2 An exception is provided by some versions of Affleck-Dine Baryogenesis [10] where a baryon excess 
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Thus, because Eq. (2) does not agree with the observed value given in Eq. (1), one is led 
to the interesting conclusion that a primordial matter-antimatter asymmetry must have 
existed at temperatures of Oil GeV) in the Universe. The observed value for %/n 7 and 
the lack of antimatter in the Universe are manifestations of this primordial asymmetry. 
Hence, in reality, the ratio t]b is, in effect, a measure of the number density of matter 
minus that of antimatter relative to the photon number density: 

= n B -n B = 6 l ± Q 3 x 1Q _ 10 

It is interesting to consider the physical origins of this primordial matter-antimatter 
asymmetry. From the seminal work of Sakharov [11] one knows that, under certain con- 
ditions which we will amplify later on in this article, this asymmetry can be generated 
by physical processes. In this review we will focus on Leptogenesis - the creation of a 
primordial lepton-antilepton asymmetry - as the root source for the observed baryon- 
antibaryon asymmetry of Eq. (3) [12]. In our view, Leptogenesis provides the most com- 
pelling scenario for generating the observed baryon asymmetry in the Universe. In par- 
ticular, because Leptogenesis is closely linked with parameters in the neutrino sector 
that can be eventually determined experimentally, this scenario can be tested and can 
be either confirmed or ruled out by data. 

The plan of this review is as follows. In Section 2 we discuss the theoretical condi- 
tions necessary for producing a primordial matter- antimatter asymmetry in the Universe 
and explain how, through a mechanism first discussed by Kuzmin, Rubakov, and Sha- 
poshnikov, [13], it is possible to turn a lepton asymmetry into a baryon asymmetry. In 
Section 3 we review existing experimental information on the neutrino sector, as well 
as the seesaw mechanism [14] that provides a theoretical framework for understanding 
this data. Section 4 discusses thermal Leptogenesis and contains the main quantitative 
results of the review. In particular, we show in this Section that the observed value for 
i]b obtained from Leptogenesis significantly constrains low energy neutrino properties, 
and vice versa. In Section 5 we turn to dark matter and discuss how supersymmetric 
candidates for dark matter are significantly constrained if thermal Leptogenesis is the 
source of the observed baryon asymmetry in the Universe. The constraint arises from 
the overproduction of gravitinos. Section 6 discusses nonthermal Leptogenesis and other 
nonthermal processes that can lead to Baryogenesis. Finally, we present our conclusion 
and summary of results in Section 7. 

is produced by the decay of a scalar field very late in the history of the Universe, which reheats the 
Universe to temperatures of the 0(100 MeV). 
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2 Theoretical Foundations 



2.1 Sakharov's Conditions for Baryogenesis 

In 1967, Sakharov [11] considered the consequences of the hypothesis that the observed 
expanding Universe originated from a superdense initial state with temperature of order 
the Planck mass, ~ Mp. Since he could not imagine how, starting from such an initial 
state, one could obtain a macroscopic separation of matter and antimatter, he concluded 
that our Universe contains today only matter. That is, the Universe evolved from an 
initial state even under charge conjugation to a state odd under charge conjugation 
today He then realized that in an expanding Universe a matter-antimatter asymmetry 
could be generated dynamically, if C, CP, and baryon and lepton number were violated, 
and these processes were out of thermal equilibrium. 

Sakharov also described a concrete model for Baryogenesis. He proposed as the origin 
for the baryon and lepton asymmetry the CP-violating decays of maximons, hypothetical 
neutral spin zero particles with mass of order the Planck mass. Their existence leads to 
a departure from thermal equilibrium already at temperatures T ~ Mp, where a small 
matter-antimatter asymmetry is then generated. An unavoidable consequence of this 
model is that protons are unstable and decay. However, the proton lifetime in Sakharov's 
model turned out to be unobservably long, r v > 10 50 years. 

During the past four decades many models of Baryogenesis have been proposed, 
demonstrating that the conditions Sakharov spelled out to allow Baryogenesis to take 
place are quite readily satisfied in the Standard Model of particle physics and its exten- 
sions. There are, however, significant differences among the various mechanisms suggested 
for producing the baryon asymmetry. Grand Unified Theories (GUTs) have been of par- 
ticular importance for the development of realistic models of Baryogenesis [15]. These 
theories provide natural heavy particle candidates, whose decays can be the source of 
the baryon asymmetry. However, in general, the simplest GUT models based on SU(5) 
lead to a creation of a (B+L)-asymmetry, with a vanishing asymmetry for B-L. As will 
be made clear below, a (B+L)-asymmetry generated at the GUT scale eventually gets 
erased by sphaleron processes. In Leptogenesis, heavy Majorana neutrinos required by 
the seesaw mechanism [14] serve to trigger Baryogenesis. Because B-L is violated, the 
erasure present in GUTs is avoided. In principle, Electroweak Baryogenesis [16] is also 
an attractive possibility, as the relevant parameters could then be tested in collider ex- 
periments. However, in general, the electroweak phase transition is not sufficiently out of 
equilibrium to generate an asymmetry of the magnitude observed in the Universe. Finally, 
in supersymmetric theories the baryon and lepton number stored in scalar expectation 
values can also lead to Baryogenesis, through the, so-called, Affleck-Dine mechanism [10], 
which will be discussed in some detail in Section 6. 
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2.2 B + L Violation in the Standard Model 



Due to the chiral nature of the electroweak interactions, baryon and lepton number are 
not conserved in the Standard Model [17]. The divergence of the B and L currents, 

= \ Yl (QLl^L+URl^UR + dji-f^dj^ , (4) 

generations 

Ji = (hluh + eRl^ , (5) 

generations 

is given by the triangle anomaly, 

= 32^ (-9 2 K^ U + 9' 2 B UU B^) . (6) 

Here Nf is the number of generations, and and B u are, respectively, the SU (2) and 
U(l) gauge fields with gauge couplings g and g'. 

As a consequence of the anomaly, the change in baryon and lepton number is related 
to the change in the topological charge of the gauge field, 

B(t f ) - B(U) = J' 1 dt J d 3 xd>*J* 

= Nf[N cs (t f ) - N cs (ti)] , (7) 

where 

N cs (t) = tJ^ / d 3 xe tjk e IJK W n W J iW Kk . (8) 

For vacuum to vacuum transitions W I% is a pure gauge configuration and the Chern- 
Simons numbers N cs {ti) and N cs (tf) are integers. 

In a non-abelian gauge theory there are infinitly many degenerate ground states, 
which differ in their value of the Chern-Simons number, AiV cs = ±1, ±2, . . .. The corre- 
ponding points in field space are separated by a potential barrier whose height is given 
by the so-called sphaleron energy E sph [18]. Because of the anomaly, jumps in the Chern- 
Simons number are associated with changes of baryon and lepton number, 

AB = AL = N f AN cs . (9) 

Obviously, in the Standard Model the smallest jump is AB = AL = ±3. 

In the semiclassical approximation, the probability of tunneling between neighboring 
vacua is determined by instanton configurations. In the Standard Model, SU (2) instan- 
tons lead to an effective 12-fermion interaction 

O b +l = II (QuQuQulu) , (10) 

i=l. ..3 
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which describes processes with AB = AL = 3, such as 



u c + d c + c c -> d + 2s + 2b + t + v e + ^ + i/ T . (11) 
The transition rate is determined by the instanton action and one finds [17] 



C 47T 



= £>(l(T 165 ) . (12) 



Because this rate is extremely small, [B + L)-violating interactions appear to be com- 
pletely negligible in the Standard Model. However, this picture changes dramatically 
when one is in a thermal bath. 



2.3 Sphalerons and the KRS Mechanism 

As emphasized in the seminal paper of Kuzmin, Rubakov and Shaposhnikov [13], in the 
thermal bath provided by the expanding Universe one can make transitions between the 
gauge vacua not by tunneling, but through thermal fluctuations over the barrier. For 
temperatures larger than the height of the barrier, the exponential suppresion in the 
rate provided by the Boltzmann factor disappears completely. Hence (B+L)-violating 
processes can occur at a significant rate and these processes can be in equilibrium in the 
expanding Universe. 

The finite-temperature transition rate in the electroweak theory is determined by 
the sphaleron configuration [18], a saddle point of the field energy of the gauge-Higgs 
system. Fluctuations around this saddle point have one negative eigenvalue, which allows 
one to extract the transition rate. The sphaleron energy is proportional to Vf(T), the 
finite-temperature expectation value of the Higgs field, and one finds 

E sph (T) ~ jV F (T) . (13) 

Taking translational and rotational zero-modes into account, one obtains for the transi- 
tion rate per unit volume in the Higgs phase [19] 

Vb + l - K M w -0E mh (T) / 1 A) 

where (3 = 1/T, M w = g 2 vp(T)/2 and k is some constant. 

Extrapolating this semiclassical formula to the high-temperature symmetric phase, 
where v F (T) = 0, and using for M w the thermal mass, M w ~ g 2 T, one expects in this 
phase Tb+l/V ~ (aT) 4 . However, detailed studies have shown that this naive extrapo- 
lation from the Higgs to the symmetric phase is not quite correct. The relevant spatial 



scale for non-perturbative fluctuations is the magnetic screening length ~ l/(g 2 T), but 
the corresponding time scale turns out to be 1 / (g 4 T In g^ 1 ), which is larger for small 
coupling [20,21]. As a consequence one obtains for the sphaleron rate in the symmetric 
phase 

Tb+l/V ~ a 5 In a _1 T 4 . (15) 

It turns out that the dynamics of low-frequency gauge fields can be described by 
a remarkably simple effective theory, derived by Bodeker [21]. The color magnetic and 
electric fields satisfy the equation of motion 

DxB = aE-(. (16) 

Here ( is Gaussian noise, a random vector field with variance 

{( i (x)( j (x')) = 2cT6 ij 6\x-x'). (17) 

These equations define a stochastic three-dimensional gauge theory. The parameter a is 
the 'color conductivity', a = m 2 D /{?> r y), where nip ~ gT is the Debye screening mass and 
7 ~ g 2 Tln(l/g) is the hard gauge boson damping rate. To leading-log accuracy one has 
1/cr ~ \iig~ 1 . A next-to-leading order analysis yields for the sphaleron rate [22] 




The overall coefficient has been determined by a numerical lattice simulation [23]. From 
Eq. (|18|) one easily obtains the temperature range where sphaleron processes are in 
thermal equilibrium: 

Tew ~ 100 GeV < T < T sph ~ 10 12 GeV . (19) 

The effective theory describing topological fluctuations of the gauge field in the high- 
temperature phase is valid for small coupling, g <C 1. Yet for T EW < T < T sph ~ 10 12 GeV 
one has g = 0(1). This implies that the electric screening lenth l/(gT) and the magnetic 
screening length 1 / (g 2 T) are not well separated and that nonperturbative corrections to 
the sphaleron rate, Eq. (|18j) . may be large. This will modify the temperature range given 
in Eq. ()19j) . but one expects that the qualitative picture of fluctuations in baryon and 
lepton number in the high-temperature phase of the Standard Model will not be affected. 

2.4 Electroweak Baryogenesis and its Experimental Con- 
straints 

An important ingredient in the theory of Baryogenesis is related to the nature of the elec- 
troweak transition from the high-temperature symmetric phase to the low-temperature 
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Higgs phase. Because in the Standard Model baryon number, C and CP are not con- 
served, it is conceivable that the cosmological baryon asymmetry could have been gener- 
ated at the electroweak phase transition [13], provided that this transition is of first-order, 
because then there is also the necessary departure from thermal equilibrium. This pos- 
sibility has stimulated much theoretical activity during the past years to determine the 
phase diagram of the electroweak theory. 

Electroweak Baryogenesis requires that the baryon asymmetry generated during the 
phase transition is not erased by sphaleron processes afterwards. 3 This leads to a con- 
dition on the jump of the Higgs vacuum expectation value Vp = V WH at the critical 
temperature [24]: 

^fP>l. (20) 

The strength of the electroweak transition has been studied by numerical and analytical 
methods as function of the Higgs boson mass. For the SU{2) gauge-Higgs model one 
finds from lattice simulations as well as perturbative calculations that the lower bound 
of Eq. (f2T)J) is violated for Higgs masses above 45 GeV [25]. 4 Because the present lower 
bound from LEP on the Higgs mass is 114 GeV [26], it is clear that the electroweak 
transition in the Standard Model is too weak for Baryogenesis. However, for special 
choices of parameters or by adding singlet fields, in certain circumstances supersymmetric 
extensions of the Standard Model have a sufficiently strong first-order phase transition 
to allow Electroweak Baryogenesis to take place [27]. 

For large Higgs masses, the nature of the electroweak transition is dominated by 
nonperturbative effects of the SU(2) gauge theory at high temperatures. At a critical 
Higgs mass m c H = 0{My/\ an intriguing phenomenon occurs: The first-order phase 
transition turns into a smooth crossover [28-30], as expected on general grounds [25]. At 
the endpoint of a critical line of first-order transitions, which is reached for itlh = m c H , 
the phase transition is of second order [31]. 

The value of the critical Higgs mass can be estimated by comparing the W-boson 
mass My/ in the Higgs phase with the magnetic mass tusm in the symmetric phase. This 
yields m c H ~ 74 GeV [32]. Numerical lattice simulations have determined the precise 
value m c H = 72.1 ± 1.4 GeV [33]. The analytic estimate of the critical Higgs mass can 
be generalized to supersymmetric extensions of the Standard Model, where one finds 
m c h < 130 . . . 150 GeV [34], which is still compatible with the present experimental lower 
bound. 



3 The produced asymmetry will be erased if, after the phase transition, (B+L)-violating processes are 
in equilibrium. 

4 For Higgs masses below 50 GeV, the Higgs model provides a good approximation for the full Standard 
Model. 
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2.5 The Relation Between Baryon and Lepton Asymmetries 



In a weakly coupled plasma, one can assign a chemical potential fi to each of the quark, 
lepton and Higgs fields. In the Standard Model, with one Higgs doublet H and Nf 
generations one then has 5Nf + 1 chemical potentials. 5 For a non-interacting gas of 
massless particles the asymmetry in the particle and antiparticle number densities is 
given by 

_ _ = gT^ \ PfH + O ((Pfii) 3 ) , fermions , 
Ul Ui Q { 2/3//; + O ((/9/ii) 3 ) , bosons . [ ' 

The following analysis is based on these relations for /3fii <C 1. However, one should keep 
in mind that the plasma of the early Universe is very different from a weakly coupled 
relativistic gas, owing to the presence of unscreened non-abelian gauge interactions, where 
nonperturbative effects are important in some cases. 

Quarks, leptons and Higgs bosons interact via Yukawa and gauge couplings and, in 
addition, via the nonperturbative sphaleron processes. In thermal equilibrium all these 
processes yield constraints between the various chemical potentials [35]. The effective 
interaction of Eq. induced by the SU(2) electroweak instantons implies 

(3/v + Wi) = • (22) 

i 

One also has to take the 377(3) Quantum Chromodynamics (QCD) instanton processes 
into account [36], which generate an effective interaction between left-handed and right- 
handed quarks. The corresponding relation between the chemical potentials reads 

(2/igi - \i u i ~ fJ>di) = . (23) 

i 

A third condition, valid at all temperatures, arises from the requirement that the total 
hypercharge of the plasma vanishes. From Eq. (|21|) and the known hypercharges one 
obtains 

Yl (^qi + 2/ini - A*di - m ~ Hei + J^V^ = • (24) 

The Yukawa interactions, supplemented by gauge interactions, yield relations be- 
tween the chemical potentials of left-handed and right-handed fermions, 

AV - HH- Hdg = , [Iqi + {J>H Huj = , Uli ~ HH ~ Vej = . (25) 

These relations hold if the corresponding interactions are in thermal equilibrium. In the 
temperature range 100 GeV < T < 10 12 GeV, which is of interest for Baryogenesis, 

5 In addition to the Higgs doublet, the two left-handed doublets qt and li and the three right-handed 
singlets m, d%, and e; of each generation each have an independent chemical potential. 
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this is the case for gauge interactions. On the other hand, Yukawa interactions are in 
equilibrium only in a more restricted temperature range that depends on the strength of 
the Yukawa couplings. In the following we shall ignore this complication which has only 
a small effect on our discussion of Leptogenesis. 

Using Eq. the baryon number density % = gBT 2 /6 and the lepton number 

densities = LigT 2 /6 can be expressed in terms of the chemical potentials: 

B = (2/i qi + /i ui + /i di ) , (26) 

i 

Li = 2/iu + fi ei , L = ^U . (27) 

i 

Consider now the case where all Yukawa interactions are in equilibrium. The asym- 
metries Li — B/Nf are then conserved and we have equilibrium between the different 
generations, jiu = /ij, n q i = n q) etc. Using also the sphaleron relation and the hyper- 
charge constraint, one can express all chemical potentials, and therefore all asymmetries, 
in terms of a single chemical potential that may be chosen to be /i;, 

2N f + 3 6N f + 1 2N f - 1 

^ ~ 6iV> + 3^' ^-'GNf + Z*' liu ~ 6^ + 3^' 

1 ANf 
N = ^=^3^1. (28) 

The corresponding baryon and lepton asymmetries are 

AN f UN 2 + §N f , 

b = L m , L = — ^ -m . 29 

3 P QNf + 3 P v ; 

This yields the important connection between the B, B — L and L asymmetries [37] 

B = c s {B-L); L = {c s -1)(B-L) , (30) 

where c s = (8Nf + A)/(22Nf + 13). The above relations hold for temperatures T 3> Vp. 
In general, the ratio BjiB — L) is a function of Vp/T [38]. 

The relations (J3*U|) between B-, (B-L)- and L-number suggest that (B-L)-violation is 
needed in order to generate a B-asymmetry. 6 Because the (B-L)-current has no anomaly, 
the value of B-L at time tf, where the Leptogenesis process is completed, determines the 
value of the baryon asymmetry today, 

B(t ) = c s {B-L){t f ) . (31) 



6 In the case of Dirac neutrinos, which have extremely small Yukawa couplings, one can construct 
Leptogenesis models where an asymmetry of lepton doublets is accompanied by an asymmetry of right- 
handed neutrinos such that the total L-number is conserved and the (B-L)-asymmetry vanishes [39] . 
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On the other hand, during the Leptogenesis process the strength of (B-L)-, and therefore 
L- violating interactions can only be weak. Otherwise, because of Eq. ()30j) . they would 
wash out any baryon asymmetry. As we shall see in the following, the interplay between 
these conflicting conditions leads to important constraints on the properties of neutrinos. 

3 Experimental and Theoretical Information on the 
Neutrino Sector 

3.1 Results from Oscillation Experiments 

The search for neutrino mass has a long history [40]. Positive results are now provided 
by neutrino oscillation experiments. The allowed values for the mass-squared differences 
Am| and the mixing angles 6^ at the 3a level for three generations of neutrinos are 
summarized below [26]: 

sin 2 2# 23 = 0.92 - 1, |Am|,| = (1.2 - 4.8) x 10~ 3 eV 2 (atmospheric v) (32) 

sin 2 2# 12 = 0.70 - 0.95, |Am 2 2 | = (5.4 - 9.5) x 10" 5 eV 2 (solar v). (33) 
The CHOOZ experiment [41] gives only an upper limit on the remaining mixing angle 

sin 2 2# 13 = - 0.17 (CHOOZ v). (34) 

The LSND experiment [42] reports neutrino oscillations from to u e . The mixing angle 
and mass-squared difference inferred from this experiment are sin 2 2# = 0.003 — 0.03 and 
| Am 2 1 = 0.2 — 2 eV 2 . This parameter region is almost excluded by negative results from 
a comparable experiment by the KARMEN collaboration [43], but there still remains a 
narrow region allowed at the 90% CL. 

It is very difficult to explain all the above data from neutrino oscillation experiments 
within the three neutrino framework. Indeed, to accommodate the LSND data one must 
introduce at least one sterile neutrino [44], or make the radical assumption that CPT 
is not conserved [45]. In this review, for simplicity we will disregard the data from the 
LSND experiment. 

3.2 Information from /3-decay, 2/3-decay and Cosmology 

The oscillation experiments discussed in the previous subsection are only sensitive to 
mass-squared differences. In this subsection we quote results from direct searches for the 
absolute values of neutrino masses. 
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The direct laboratory limits on the neutrino masses are summarized as follows [26] : 

m„ e < 2.5 eV ; (35) 

m„ M < 170 keV ; (36) 

m Vr < 18 MeV . (37) 

The v e mass measurements use the decay of tritium, 3 H — > 3 He + v e + e~, which has 
a small Q value, Q = 18.6 keV, and looks at the electron spectrum near the end point 
in the Kurie plot. The limit on the mass is obtained from the two-body kinematics 
of the pion decay, 7r + — > /i + + v^. Finally, the limit on the v T mass is obtained from 
measurements of the invariant mass distribution of 3n and 5n systems in the r decays, 

T — > 3(5)7T + V T . 

Neutrinoless double j3 decay experiments provide a bound on an element of the 
Majorana mass matrix, m Uee [46]. The best limit comes from the 76 Ge results. Because 
the calculation of the double (3 decay rate is model dependent, we quote a range for this 
bound: 

m^ e < 0.3 - 0.8 eV . (38) 

One can derive a stringent upper limit on the sum of neutrino masses from cosmology. 
In the early Universe neutrinos were in thermal equilibrium with radiation and one can 
infer their number density today to be n« ~ 110 cm -3 for each neutrino species. Although 
the contribution of massless neutrinos to the Universe's energy density is negligible, 
the contribution of massive neutrinos could be important. By requiring that the energy 
density of massive neutrinos does not exceed that of dark matter, one obtains the bound: 

Y / m i <30h 2 eV, (39) 

i 

where h is the Hubble constant in units of 100 km s^Mpc -1 [h = 0.7ll^| [4]]. 

Even if neutrinos satisfy the above limit, massive neutrinos would affect the formation 
of cosmic structure, because the free streaming of neutrinos suppresses density fluctua- 
tions at small scales. The normalization of large- and small-scale fluctuations constrains 
the contribution from neutrinos. Recent detailed analyses [47] lead to the bound: 

^mi<0.65eV. (40) 

i 

It is interesting that a similar constraint, J2i m i < 2.0 eV, has been obtained by using 
the cosmic microwave background data alone [48]. 
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3.3 The Seesaw Mechanism 



If there are right-handed neutrinos, then neutrinos can have a Dirac mass much as quarks 
and charged leptons do. However, if this is the only source for their mass, it is not easy 
to find a natural reason for the very small mass of neutrinos. With only a Dirac mass 
term, the smallness of neutrino masses needs to be ascribed to having very tiny Yukawa 
coupling constants h (h ~ 1CT 13 gives a neutrino mass m u ~ 0.01 eV). Although it is 
possible to imagine mechanisms that result in very small Dirac masses for neutrinos, 7 the 
seesaw mechanism, which entails Majorana masses, provides a natural explanation for 
the smallness of neutrino masses in theories of unification at ultra-high energy scales [14]. 

To appreciate this point, we should first note that the Standard Model does not 
require neutrinos to be massive, because right-handed neutrinos are not necessary for 
the electroweak theory. Without right-handed neutrinos, these particles may acquire 
their mass only from what are called irrelevant operators - operators such as UHH 
with dimension greater than four. These operators can give rise to a Majorana mass for 
neutrinos in theories with a cutoff. However, in the limit of an infinite cutoff, neutrino 
masses vanish. It should be noted that if one adopts the Planck scale as the Wilson 
cutoff for the Standard Model, one finds neutrino masses to be at most 10~ 5 eV. Thus, 
the observed mass \J Aml tm ~ 0.05 eV is unable to be explained within the Standard 
Model. 

If one considers possible extensions of the Standard Model gauge group, it is natural 
to consider a gauge group G whose rank is at least 5, since the rank of the Standard 
Model group, SU (3) x SU (2) x U(l)y , is 4. This new group G may contain an extra U(l) 
as a subgroup, in addition to the Standard Model gauge group. The simplest candidate 
for the extra £7(1) is a B-L gauge symmetry, because we know that the global U(1) B -l 
does not have any anomaly due to the Standard Model gauge interactions. However, 
when one gauges this B-L symmetry, the theory does have a self-anomaly of the B-L 
interactions. That is, the triangle anomaly of [U(1)b-l] 3 is non-vanishing. A crucial 
point, however, is that this anomaly is cancelled by introducing right-handed neutrinos. 
This also cancels the mixed gravitational/B-L anomaly. Thus, right-handed neutrinos 
are required for consistency of the theory! A famous example which includes U(1)b~l is 
provided by S'O(IO) grand unification. But our argument is more general. For instance, 
the string brane world predicts many £/(l)'s and it is quite natural to consider some 
of them to be anomaly free and survive as low-energy (compared with the string scale) 
gauge symmetries. 

It is usually assumed that the unification group G is broken down to the Standard 

7 At the end of this subsection we outline a recent approach that may explain how such extremely 
small Yukawa coupling constants might arise in a higher dimensional theory. 
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Model group at high energies. Then, the right-handed neutrinos naturally obtain large 
Majorana masses, because they are singlets of the Standard Model and there is no unbro- 
ken symmetry to protect them from acquiring a large Majorana mass. In this article we 
shall denote heavy Majorana neutrinos as N. Then, the masses of the neutrinos written 
as a matrix take a simple form: 

m 
m T M 



(41) 



Here m is the Dirac mass matrix between the left-handed neutrinos v and the heavy 
Majorana neutrinos N, which is of order the electroweak scale, while M is the Majorana 
mass matrix of the heavy neutrinos N. For three generations, both M and m are 3x3 
matrices. Integrating out the heavy Majorana neutrinos N leads to a small neutrino mass 
via the seesaw mechanism [14]. For one neutrino generation one simply has that: 

2 

m v ~ — . (42) 

We see from the above that a small neutrino mass is a reflection of the ultra-heavy 
mass of the heavy neutrino N. The observed neutrino mass \J Am^ ~ 0.05 eV implies 
M ~ 10 15 GeV, which is very close to the Grand Unification scale. Thus, effectively, the 
small neutrino masses provide a window to new physics at an ultra-high energy scale. 

One may question why the unification group G, or the U(1)b-l symmetry, should 
be broken at such a high energy scale. Perhaps the answer to this question, as we shall 
amplify in this review, is because, otherwise, the baryon number in the Universe would 
be too small for us to exist! It turns out that if the Universe's baryon number were to 
be two orders of magnitude below the present observed value galaxies would not be 
formed [49]. One of the purposes of this review article is to explain in some detail this 
fundamental point. 



3.3.1 Small Neutrino Yukawa Couplings from Higher Dimensional Theories 

To explain how one can generate small Dirac masses for neutrinos, consider a theory 
described in (4+l)-dimensional space-time, while our world is on a (3+l)-dimensional 
hyperplane - a so-called D3 brane. The Einstein action of gravity in five-dimensional 
space-time is given by 

where M* is the gravitational scale in five-dimensional space-time, and g$ and IZ5 are 
the metric and the scalar curvature, respectively. We assume that the fifth dimension is 
compactified to a space of radius L, and consider the metric to be 



d 2 s = g^dx^dx" — dy 2 , 
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(44) 



where g^ v is the metric in four-dimensional space-time. The integration over dy leads to 
the four-dimensional action 

Then, the Planck scale, M P ~ 1.2 x 10 19 GeV, in four-dimensional space-time is given 
by 

Mp = MlL. (46) 

One can get the observed value for M P even for M* = 1 TeV by taking a very large L. 
The weakness of gravity in these theories is the result of having a large compactification 
scale in the fifth dimension [50]. 

Let us now assume that all the Standard Model particles reside on a D3 brane at 
the boundary y = and that the right-handed neutrino lives in the five-dimensional 
bulk [51]. Then the action involving the right-handed neutrinos is given by 

S = M*J d 4 x £ dyy/=f 5 N R i$N R + J d 4 x £ d yv ^g~ 5 hN R £ L H5(y) + h.c. , (47) 

Here, £l and H are SU(2)^ doublets of the left-handed leptons and of the Higgs boson, 
respectively. One obtains the action in four dimensions by integrating over dy, and one 
finds: 

S 4 = M*lJ d A x^f A N R i$N R + J d A x^/^lhN R l L H + h.c. . (48) 

After renormalizing the wave function of N R so that it has a canonical kinetic term, one 
finds that the Yukawa coupling constant is suppressed by l/y/M*L = M*/M P . Thus in 
this model one obtains an effective Yukawa coupling constant h e g ~ 10~ 13 (corresponding 
to a neutrino Dirac mass m v ~ 0.01 eV ) for h — 1 and M* ~ 10 3 TeV. 

This model, however, has a serious drawback. There is no symmetry to protect the 
right-handed neutrinos from acquiring a Majorana mass, because they are singlets of the 
Standard Model gauge group. A solution to this problem may be found by imposing a 
U(1)b-l gauge symmetry in the five-dimensional bulk. As long as the B-L gauge symme- 
try is exact, the right-handed neutrinos cannot have a Majorana mass because this mass 
carries a no n- vanishing B-L charge. If this symmetry is exact, the corresponding gauge 
boson is completely massless. However, this may not cause any phenomenological diffi- 
culties at low-energies, because the B-L gauge coupling constant must also be extremely 
suppressed. 8 But, as we explained in Section 2, if the B-L gauge symmetry is exact, it is 
very difficult to account for the baryon-number asymmetry in the present Universe. 

8 Thc B-L gauge coupling constant as-L is constrained to be olb-l < 10~ 21 a O m- This bound comes 
from the empirical limits of the electromagnetic charges for the neutron and the neutrino. That is, 
Q n = (-0.4 ± 1.1) x 1(T 21 [52] and Q v = (0.5 ± 2.9) x 10~ 21 [53]. The present model suggests a B -L ^ 
ijj^) 2 x a cm ~ 10~ 25 x a om , which may be in an interesting region for future experiments. 
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3.4 CP-Violating Phases at Low and High Energies in the 
Lepton Sector 



In the Standard Model, the Lagrangian for the lepton sector, augmented by including 
right handed neutrinos, is given by 

C = lui^Li + e Ri i$e Ri + N Ri i$N Ri 

■ .(;;,< H.h.jU" + h i:j N Ri £ Lj H - )^M tJ N Rl N R] + h.c. , (49) 

where i,j = {1 — 3} are the family-number indices. We adopt, without losing generality, 
a basis where the matrices fa and M^- are diagonal. The Yukawa matrix hij in this basis 
is in general complex and thus has CP-violating phases. Because for three families, the 
matrix has 9 complex parameters, we have 9 possible CP-violating phases. However, 
three of these phases can be absorbed into the wave function of £ L and hence 6 CP- 
violating phases remain physically relevant. These are known as high-energy phases, 
because they enter in the full theory. 9 

Let us now discuss the CP-violating phases at low energies. To do that we first need 
to integrate out the heavy Majorana neutrinos, iVj. Doing so the effective Lagrangian for 
the lepton sector reduces to: 

1 - H 2 

£ e ff = iLitfliA + e Ri i^e Ri + fae Ri l Li H ] + - V hJ k h kj £ Li £ Lj — + h.c. . (50) 

2 k M k 

The last term can be rewritten as 

1 H 2 

- ^m Uij £ Li £ Lj -^, (51) 

so that all low-energy phases appear in the mass matrix of light neutrinos. Because the 
neutrino mass matrix is symmetric, it has 6 complex parameters, and hence one has 6 
possible CP-violating phases. However, as before, 3 of these 6 phases can be absorbed 
into the wave function of £l- Therefore, there remain only three physical low energy CP- 
violating phases [54] . Because they are different in number, it is unfortunately very diffi- 
cult to establish a direct link between the low-energy and the high-energy CP-violating 
phases [55]. 

Furthermore, in practice, it is not possible to measure all three low-energy phases. 

One of these three phases can be measured by neutrino oscillation experiments, while 

neutrinoless double f3 decay, if it were to be observed, would provide information on 

another phase. However, the remaining phase is undetermined. In other words, one cannot 

9 In particular, as we will show in the next Section, the phase contributing to the generation of the lep- 
ton asymmetry in the decay of Ni is a combination of these high-energy phases, given by yj lm[(hh^)n] 2 . 
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perform a complete experiment to determine the neutrino mass matrix. Nevertheless, if 
the neutrino mass matrix were to have an extra constraint, one may be able to determine 
all matrix elements oim v . This constraint must be independent of the frame of the family 
basis. One example of such a constraint is the requirement that det(m„) = 0. In this 
case, we have only 7 independent physical parameters including the phases in the neutrino 
mass matrix [56], which can be determined in principle in future experiments. 10 



4 Thermal Leptogenesis 

4.1 Lepton Number Violation and Leptogenesis 

As we discussed above, lepton number violation is most simply realized by adding right- 
handed neutrinos to the Standard Model. Their existence is predicted by all extensions 
of the Standard Model containing B-L as a local symmetry and allows for an elegant 
explanation of the smallness of the light neutrino masses via the seesaw mechanism [14]. 

The most general Lagrangian for couplings and masses of charged leptons and neu- 
trinos is given in Eq. ([49)1 . The vacuum expectation value of the Higgs field, (H) = vp, 
generates Dirac masses m e and vtld for charged leptons and neutrinos, m e = fvp and 
nip = hvF, which are assumed to be much smaller than the Majorana masses M. This 
yields the light and heavy neutrino mass eigenstates 

u~V?u L + viV; , N~N R + N C R , (52) 

with masses ^ 

m v ~ -Vjm T D —m D V v , m N ~ M . (53) 

In a basis where the charged lepton mass matrix m e and the Majorana mass matrix M 
are diagonal, V v is the mixing matrix in the leptonic charged current. 

The right-handed neutrinos can efficiently erase any pre-existing lepton asymmetry 
at temperatures T > M, but they can also generate a lepton asymmetry by means of 
their out-of-equilibrium decays at temperatures T < M. This asymmetry is then partially 
transformed into a baryon asymmetry by sphaleron processes. This is the Leptogenesis 
mechanism proposed by Fukugita and Yanagida [12]. 

The decay width of the heavy neutrino iVj at tree level reads, 

V Di = T(N, t ^ H + £ L ) + T (Ni ^H^ + i{) = -^-(hh^M . (54) 

v 7 S7T 

10 As will be shown in Section 6, Affleck-Dine Leptogenesis suggests a constraint, m Vl ~ 10 -9 eV and 
hence det(m„) ~ 0. 
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Figure 1: Tree level and one- loop diagrams contributing to heavy neutrino decays whose 
interference leads to Leptogenesis. 

Once the temperature of the universe drops below the mass Mi, the heavy neutrinos are 
not able to follow the rapid change of the equilibrium distribution. Hence, the necessary 
deviation from thermal equilibrium ensues as a result of having a too large number 
density of heavy neutrinos, compared to the equilibrium density. Eventually, however, the 
heavy neutrinos decay, and a lepton asymmetry is generated owing to the presence of CP- 
violating processes. The CP asymmetry involves the interference between the tree-level 
amplitude and the one-loop vertex and self-energy contributions (see Fig. (1)). In a basis, 
where the right-handed neutrino mass matrix M is diagonal, one obtains [57] for the CP 
asymmetry parameter E\ assuming hierarchical heavy neutrino masses (M\ <C M2,Ms): 



£1 ~ Im 

" 167T (/ i / i t) n 2- 



hh^ 2 

' il 



In the case of mass differences of order the decay widths, one obtains a significant en- 
hancement from the self-energy contribution [58], although the influence of the thermal 
bath on this effect is presently unclear. 

The CP asymmetry of Eq. (|55jl can be obtained in a very simple way by first inte- 
grating out the heavier neutrinos N 2 and N 3 in the leptonic Lagrangian. This yields 

£ e J f = h x ^T Rx t L3 H - l -Mj^ Rl N RX + ^ Vij £ u He Lj H + h.c. , (56) 

with 

3 1 

Vij = h Ikiri- h kj ■ (57) 

k=2 M k 

The asymmetry e\ is then obtained from the interference of the Born graph and the one- 
loop graph involving the cubic and the quartic couplings. This includes automatically 
both, vertex and self-energy corrections [59] and yields an expression for e\ directly in 
terms of the light neutrino mass matrix: 



3 Mi 

^2 



167r {hW)nv F 

The CP asymmetry then leads to a (B-L)-asymmetry [12] 



Im (h*m u hJ) . {hi 



Y B _ L c -Y L = = _„£l . (59) 

s g* 
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Here s is the entropy and, in the present epoch s = 7.04n 7 , whereas g* ~ 100 is the 
number of degrees of freedom in the plasma. The factor k < 1 in the above takes into 
account the effect of washout processes. As we shall discuss below, in order to determine 
k one has to solve the Boltzmann equations. 

Early studies of Leptogenesis were partly motivated by trying to find alternatives 
to Electroweak Baryogenesis, which did not seem to produce a big enough asymmetry. 
Some extensions of the Standard Model were considered and, in particular, in the simple 
case of hierarchical heavy neutrino masses the observed value of the baryon asymmetry 
is naturally obtained with B-L broken at the unification scale, M GUT ~ 10 15 GeV. The 
corresponding light neutrino masses are then very small, < m 3 ~ 0.1 eV, and the 
typical parameters for the necessary CP asymmetry and the Baryogenesis temperature 
are e 1 ~ 10~ 6 and T B ~ M 1 ~ 10 10 GeV, respectively [60, 61]. 11 Subsequently, researchers 
realized that such small neutrino masses are consistent with the small mass differences 
inferred from the solar and atmospheric neutrino oscillations. This fact has given rise 
to a strong interest in Leptogenesis in recent years, and a large number of interesting 
models have been suggested [63]. 



4.2 Departure from Thermal Equilibrium 

Leptogenesis takes place at temperatures T ~ Mi. For a decay width small compared to 
the Hubble parameter, Ti(T) < H(T), heavy neutrinos are out of thermal equilibrium, 
otherwise they are in thermal equilibrium [64]. The borderline between the two regimes 
is given by Ti = H\t=m 1 , which is equivalent to the condition that the effective neutrino 



mass 



mi = ^fhi (60) 

Mi v ' 



is equal to the 'equilibrium neutrino mass' 

167T 5 / 2 1/oW | 



Here we have used the Hubble parameter H(T) ~ 1.66 g* T 2 /Mp where g* = qsm = 
106.75 is the total number of degrees of freedom and Mp = 1.22 x 10 19 GeV is the Planck 
mass. 

It is quite remarkable that the equilibrium neutrino mass m* is close to the neu- 
trino masses suggested by neutrino oscillations, \J Am s 2 ol ~ 8 x 10~ 3 eV and y Aml tm ~ 
5 x 10~ 2 eV. This encourages one to think that it may be possible to understand the 



11 For early work based on 5O(10), see [62]. 
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Figure 2: The evolution of the Ni abundance and the B — L asymmetry for a typical 
choice of parameters, M\ = 10 10 GeV, E\ = 10~ 6 , m\ = 10 _3 eV and m = 0.05 eV. 
From [67]. 



cosmological baryon asymmetry via Leptogenesis as a process close to thermal equilib- 
rium. Ideally, AL = 1 and AL = 2 processes should be strong enough at temperatures 
above Mi to keep the heavy neutrinos in thermal equilibrium and weak enough to allow 
the generation of an asymmetry at temperatures below M\. 

In general, the generated baryon asymmetry is the result of a competition between 
production processes and washout processes that tend to erase any generated asymmetry. 
Unless the heavy Majorana neutrinos are partially degenerate, — Mi <C Mi, the 
dominant processes are decays and inverse decays of N\ and the usual off-shell AL = 1 
and AL = 2 scatterings [65,66]. 

The Boltzmann equations for Leptogenesis are 12 
dN Nl 



dz 
dN B _ L 



-(D + S) (N Nl - N%) , (62) 
f/ ^DIN^-N^-WNb.l, (63) 

where z = M\jT. The number density iVjvi and the amount of B — L asymmetry, 
Nb-L) are calculated in a portion of comoving volume that contains one photon at the 
onset of Leptogenesis, so that the relativistic equilibrium N\ number density is given by 
N^(z < 1) = 3/4. Alternatively, one may normalize the number density to the entropy 
density s and consider Yx = n-x/s. If entropy is conserved, both normalizations are 
related by a constant. 



12 We use the conventions of [67]. We have also summed over the three lepton flavours neglecting the 
dependence on the lepton Yukawa couplings [68] . 
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There are four classes of processes that contribute to the different terms in the above 
equations: decays, inverse decays, AL = 1 scatterings and AL = 2 processes mediated 
by heavy neutrinos. The first three processes all modify the Ni abundance and try to 
push it towards its equilibrium value Nff. Denoting by H the Hubble expansion rate, 
the term D = Td/(H z) accounts for decays and inverse decays, whereas the scattering 
term S = Ts/{H z) represents the AL = 1 scatterings. Decays also yield the source term 
for the generation of the B — L asymmetry, the first term in Eq. (J63j) . whereas all other 
processes contribute to the total washout term W = Tw/(H z) which competes with 
the decay source term. The dynamical generation of the Aq abundance and the B — L 
asymmetry is shown in Fig. (J2J) for typical parameters. 



4.3 Decays and Inverse Decays 

It is very instructive to consider first a simplified picture in which decays and inverse 
decays are the only processes that are effective. 13 For consistency, in this approximation 
the real intermediate state contribution to the 2 — > 2 processes has to be included. In 
the kinetic equations (Jo^j) and (|p^j) one then has to replace D + S by D and W by Wjd, 
respectively, where Wm is the contribution of inverse decays to the washout term. The 
solution for N B _ L in this case is the sum of two terms [64], 

N B -l{z) =N i B _ L e-^ dz ' WlD(z ' ) -le 1 K(z-m 1 ) . (64) 

Here the first term accounts for an initial asymmetry which is partly reduced by washout, 
and the second term describes B — L production from N\ decays. It is expressed in terms 
of the efficiency factor k [68] which does not depend on the CP asymmetry Ei, 

k{z) = y* dz' D (N Nl - N%) e -I> z " w ^ . (65) 

As we shall see, decays and inverse decays are sufficient to describe qualitatively many 
properties of the full problem. 

We will first study in detail the regimes of weak and strong washout. If just decays 
and inverse decays are taken into account, these regimes correspond, respectively, to the 
limits K 1 and K 3> 1 of the decay parameter 

^= rD( ; = o t ) = ^, (ee) 

H(z = 1) m* ' V ; 

introduced in the context of ordinary GUT baryogenesis [64]. Based on the insight into 
the dynamics of the non-equilibrium process gained from these limiting cases one can 



13 



This section follows closely [69] . 



23 



then obtain analytic interpolation formulas that describe rather accurately the entire 
parameter range. 

To proceed, let us first recall some basic definitions and formulas. The decay rate is 
given by the formula [70], 

T D (z) = T m (l\ , (67) 

where the thermally averaged dilation factor is given by the ratio of the modified Bessel 
functions K\ and K 2 , 

W K 2 {z) ' 
For the decay term D, one then obtains 



(68) 



D(z) =Kz /^J . (69) 

The inverse decay rate is related to the decay rate by 

T ID (z) = T D (z)^^-, (70) 

where iV z eq is the equilibrium density of lepton doublets. Because the number of degrees 
of freedom for heavy Majorana neutrinos and lepton doublets is the same, qn 1 — gi = 2, 
one has 

N%(z) = lz 2 K 2 (z), NP = 3 -. (71) 
This yields for the contribution of inverse decays to the washout term W: 

W ID (z) = 1 -D(z)^l. (72) 

All relevant quantities are given in terms of the Bessel functions K\ and K 2 , which can 
be approximated by simple analytical expressions. 

In the regime far out of equilibrium, K <C 1, decays occur at very small temperatures, 
z 3> 1, and the produced (B — L)-asymmetry is not reduced by washout effects. In this 
case, using Eq. (|6^|) with S = 0, the integral for the efficiency factor given in Eq. (j63|) 
becomes simply, 

K(z)~~(N^-N Nl (z)) . (73) 

The final value of the efficiency factor k,{ = k(oo) is proportional to the initial Ni 
abundance. If N[ = N^ 1 - = 3/4, then «f = 1. But if the initial abundance is zero, then 
K{ = as well. Therefore, in this region there is the well known problem that one has to 
invoke some external mechanism to produce the initial abundance of neutrinos. Moreover, 
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an initial (B-L)-asymmetry is not washed out. Thus in the regime K <C 1 the results 
strongly depend on the initial conditions and there is little predictivity. 

In order to obtain the efficiency factor in the case of vanishing initial N\- abundance, 
iV/Vi("2i) = Njf ~ 0, one has to calculate how heavy neutrinos are dynamically produced 
by inverse decays. This requires solving the kinetic equation Eq. (J62|) with the initial 
condition N r N = 0. 

Let us define a value z eq by the condition 

N Nl (z cq ) = N^(z eq ) . (74) 

Then Eq. (|62|) implies that the number density reaches its maximum at z = z eq . For 
z > z eq the efficiency factor is always the sum of two contributions, 

Kt(z) = k~(z) + k + (z) . (75) 

Here K~(z) and k + (z) correspond to the integration domains [z- x , z eq ) and [z eq , z), respec- 
tively. 

Consider first the case of weak washout, K -C 1, which implies z eq 3> 1. One then 
finds, 

N Nl ( Zeq )~^K. (76) 

It turns out that to first order in K, there is a cancellation between k + and k~ , yielding 
for the final efficiency factor 

K f (K)^^K*. (77) 

Note, that Eq. (|77j) does not hold for K > 1, because in this case z eq becomes small, and 
washout effects change the result. 

In the case of strong washout, K ^> 1, we can neglect the negative contribution k~ , 
because the asymmetry generated at high temperatures is efficiently washed out. Now 
the neutrino abundance tracks closely the equilibrium behavior. Because D oc K, one 
can solve Eq. (|62jl systematically in powers of 1/K, which yields 

D (N Nl (z) - N%{zj) = ^-Wid(z) + O (1) , (78) 

where we have used properties of the Bessel functions. From Eqs. (|63|) and (|7Hj) one 
obtains for the efficiency factor 14 

«(*) = £ f dz> ±W ID (z') e-!>"^") . (79) 
A J z cq z 
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Because k does not contribute we can take the lower limit below as z\. 
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Figure 3: Final efficiency factor when the washout term AW is neglected. From [69]. 

The integral is dominated by the contribution from a region around the value zb where 
the integrand has a maximum, which is determined by the condition 




W ID (z B ) = - (z B ) . (80) 



For K ^> 1 one has z B ^> 1, and the condition (jHUjl becomes approximately Wjd(zb) — 1, 
with Wid(z) > 1 for z < zb and Wid(z) < 1 for z > Zb- This means that the asymmetry 
produced for z < zb is essentially erased, whereas for z > z B , washout is negligible. 
Hence, the expression of Eq. (|79j) is a good approximation for the final efficiency factor. 

One finds that a rather accurate expression for zb{K) is given by 



z B {K) 




In 



\ 1024 



1 5^ 



(81) 



The integral of Eq. (|79jh which gives the final efficiency factor in terms of zb{K), is well 
approximated by 

K,{K) ~ T^kjK i 1 ~ e ~ l "" mK ) ■ (^) 
Both equations can also be extrapolated into the regime of weak washout, K <ti 1, where 
one obtains k,{ = 1 corresponding to thermal initial abundance, N* N = Nff = 3/4. At 
K ~ 3 a rapid transition takes place from strong to weak washout. Even here analytical 
and numerical results agree within 30%. For the case of zero initial N\ abundance one 
obtains an interpolation formula Kf{K) analogous to Eq. 
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The above discussion of decays and inverse decays can be extended to include AL = 1 
and AL = 2 scattering and washout processes. In the weak washout regime, K <C 1, the 
main effect is that the efficiency factor of Eq. ([77]) is enhanced to «f oc K . Relevant effects 
include scattering processes involving gauge bosons [71,72] and thermal corrections to 
the decay and scattering rates [72,73]. The range of different results is represented in 
Fig. (J3J) by the hatched region. An additional uncertainty in the weak washout regime 
is due to the dependence of the final results on the initial N\ abundance and a possible 
initial asymmetry created before the onset of Leptogenesis. 

The situation is very different in the strong washout regime. Here the final efficiency 
factor is not sensitive to the neutrino production because a thermal neutrino distribution 
is always reached at high temperatures. For mi > ~ 10 -3 eV, the effect of AL = 1 
processes on the washout is not larger than about 50%, as indicated by the hatched 
region in Fig. (|3*|) . Within these uncertainties, the final efficiency factor is given by the 
simple power law: 



Both the scale of solar neutrino oscillations, m sol = wAm 2 ol ~ 8 x 10 3 eV, and the 



scale of atmospheric neutrino oscillations, m atm = yAm 2 tm ~ 0.05 eV, are larger than 
the equilibrium neutrino mass m*. Hence, the range of neutrino masses, and therefore 
mi, indicated by neutrino oscillations lies entirely in the strong washout regime where 
theoretical uncertainties are small and the efficiency factor is still large enough to allow 
for successful Leptogenesis. 

4.4 Bounds on Neutrino Masses 

The AL = 2 processes with heavy neutrino exchange generate a contribution to the 
washout rate that depends on the absolute neutrino mass scale m 2 = m 2 + m 2 + m 2 , 



As long as AW can be neglected, the efficiency factor is independent of M\. With in- 
creasing m however, the washout rate AW becomes important and eventually prevents 
successful Leptogenesis. This leads to the upper bound on the absolute neutrino mass 
scale. [67,74] 

One can also obtain a lower bound on the heavy neutrino masses [75], because the 
CP asymmetry e\ satisfies an upper bound [75-78], which is a function of M\, rh\ and m. 
Since the rates entering the Boltzmann equations depend on the same quantities, there 
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exists for arbitrary neutrino mass matrices a maximal baryon asymmetry 7/g ax , 

Vb < vT x (^i,M u m) 

~ 0.01 £^(TOi,Ma,m) ^(to^Mito 2 ) . (85) 

Requiring the maximal baryon asymmetry to be larger than the observed one, 

rf^(m u M x ,m)>VB MB , (86) 

then yields a constraint on the neutrino mass parameters mi, Mi and m. For each 
value of m one obtains a domain in the (mi-Mi)-plane, which is allowed by successful 
baryogenesis. For m > 0.20 eV this domain shrinks to zero. One can easily translate this 
bound into upper limits on the individual neutrino masses. In a similar way, one finds a 
lower bound on Mi, the smallest mass of the heavy Majorana neutrinos. The resulting 
upper and lower bounds are [77] 

mi < 0.1 eV , Mi > 4 x 10 8 GeV , (87) 

where we have assumed thermal initial Ni abundance. The upper bound on the light 
neutrino masses holds for a normal as well as an inverted hierarchy of masses. For zero 
initial Ni abundance one obtains the more restrictive lower bound Mi > 2 x 10 9 GeV. 
For mi > m*, the baryon asymmetry is generated at the temperature Tg ~ Mi/zb < 
Mi. Hence the lower bound on the reheating temperature Tj is less restrictive than the 
lower bound on Mi. The results of a detailed analytical and numerical calculation are 
summarized in Fig. (@J). For the lower bound on the reheating temperature one finds 
Ti > 2 x 10 9 GeV [69, 72]. 15 

What is the theoretical error on the upper bound for the light neutrino masses? In 
order to answer this question one needs a full quantum mechanical treatment of Leptogen- 
esis, a challenging problem! A possible starting point is the Kadanoff-Baym equations for 
which a systematic expansion around the Boltzmann equations can be constructed [59]. 
One then has to calculate relativistic corrections, off-shell effects, 'memory effects', higher 
order loop corrections, etc. One important effect is the running of neutrino masses be- 
tween the Fermi scale and the energy scale of Leptogenesis [68,81]. Also relevant are the 
AL = 1 scattering processes involving gauge bosons [71,72]. Conceptually interesting are 
thermal corrections at large temperatures, T > Mi, which correspond to loop corrections 
involving gauge bosons and the top quark [72]. Their effect is large if thermal masses 
are treated as kinematical masses in the evaluation of scattering matrix elements. At 

15 In the supersymmetric case the CP asymmetry is enhanced but also the washout processes are 
stronger. These two effects partly compensate each other [79], leading to the slightly less stringent 
bound Ti > 1.5 x 10 9 GeV [80]. 
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Figure 4: Analytical lower bounds on Mi (circles) and 7] (dotted line) for ni\ = 0, 
^cmb _ q x ]_q-io anc j matm — 0.05 eV. The analytical results for Mi are compared 
with the numerical ones (solid lines). Upper and lower curves correspond to zero and 
thermal initial N\ abundance, respectively. The vertical dashed lines indicate the range 
(m sol ,m atm ). The gray triangle at large Mi and large rh\ is excluded by theoretical con- 
sistency. From [69]. 

sufficiently high temperatures the process N\ — > Hit is then kinematically forbidden 
whereas the process H — > N\V L is allowed by 'phase space'. On the contrary, thermal 
correction are small if they are only included as propagator effects [73]. It is important 
to clarify this issue for the treatment of non-equilibrium processes at high temperatures. 

The analysis [72] leads to the upper bound on the light neutrino masses m; < 0.15 eV. 
In [69] an upper bound of 0.12 eV has been obtained. About 0.02 eV of this difference is 
due to the different treatment of radiative corrections [81], the remaining 0.01 eV reflects 
differences in the treatment of thermal corrections. This discrepancy has to be compared 
with an uncertainty of about —0.02 eV due to 'spectator processes' [82], which have not 
been taken into account in both analyses. Hence, within the minimal seesaw model and 
the present status of theoretical calculations, the upper bound on the light Majorana 
neutrino masses is now known rather precisely. 

The main result of this section is summarized in Fig. (J3J). For mi > m*, the effi- 
ciency factor, and therefore the baryon asymmetry t)b, is independent of the initial N\ 
abundance. Furthermore, the final baryon asymmetry does not depend on the value of 
an initial baryon asymmetry generated by some other mechanism [77]. Hence, the value 
of t]b is entirely determined by neutrino properties. In this way Leptogenesis singles out 
the neutrino mass range 

10~ 3 eV < rrii < 0.1 eV . (88) 
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The firm predictions of thermal Leptogenesis open a window into the physics of the early 
universe at temperatures T B = (9(10 10 GeV), and we can ask what the implications are 
for dark matter, cosmology and particle physics. 

4.5 Triplet Models and Resonant Leptogenesis 

Measurements in neutrino physics determine the parameters of the neutrino mass matrix, 

m v = -m T D ^rn D + m^ et , (89) 

which in general contains a contribution from SU (2) triplet fields [83] in addition to the 
seesaw term generated by SU(2) singlet heavy Majorana neutrinos. So far, we have only 
considered the minimal case, m* nplct = 0. Clearly, a dominant triplet contribution would 
destroy the connection between Leptogenesis and low energy neutrino physics. 

The discovery of quasi-degenerate neutrinos with masses above the bound 0.1 eV 
would require significant modifications of minimal Leptogenesis and/or the seesaw mech- 
anism. In this case SU (2) triplet contributions to neutrino masses could be a possible 
way out [78,84,85]. Clearly, one then has no upper bound on the light neutrino masses 
anymore. Yet Leptogenesis with right-handed neutrino decays can still work yielding a 
slightly relaxed lower bound on the heavy neutrino masses. For instance, one may have 
rrii ~ 0.35 eV with M 1 > A x 10 s GeV [85]. 

Another way to reconcile quasi-degenerate light neutrinos with Leptogenesis makes 
use of the enhancement of the CP asymmetry in case of quasi-degenerate heavy neutrinos 
[86]. For instance, to raise the upper bound from 0.1 eV to 0.4 eV, a degeneracy of AM/M 
for the heavy neutrinos in the range 0.4 — 10~ 3 is required, depending on assumptions 
about the neutrino mass matrices [77, 78] . In the extreme case of 'resonant Leptogenesis' 
[71], CP asymmetries e = 0(1) are reached for degeneracies AM/M = O(10~ 10 ). In 
this case the right-handed neutrino masses may be as small as 1 TeV, which may lead to 
observable signatures at colliders. A number of models of this type have been constructed 
[87], some of which make use of the relative smallness of soft supersymmetry breaking 
terms [88]. 

4.6 Is the CP Violation in Leptogenesis Connected with the 
Low Energy CP Violation in the Neutrino Sector? 

As was shown in Section 3, the seesaw model has 6 CP-violating phases in the Yukawa 
matrix hij. Leptogenesis depends on one combination of these 6 phases. However, there 
are only 3 CP-violating phases at low energies. Hence it is impossible to determine all 
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6 phases in the theory, even if one were to measure all 3 low-energy phases. Futher- 
more, as we discussed earlier, one of these low-energy phases remains undetermined by 
experiments feasible at low energies. 

Nevertheless, the effective number of high energy CP-violating phases is reduced if 
one of the superheavy Majorana neutrinos iVj is extremely heavy and decouples from 
the seesaw system. In this case, the Yukawa couplings effectively are given by a 
2x3 matrix that contains 6 complex parameters and hence 6 phases. Three of the 6 
phases can be absorbed into the wave functions of £l and thus one is left with only 3 
CP-violating phases at high energies. In this case, the 3 low-energy CP-violating phases 
that appear in the neutrino mass matrix m v are reduced to only two physical phases, 
because det(m„) ~ 0. Although these 2 low-energy phases are, in principle, measurable 
in future experiments, this is still not enough to determine all three phases in the full 
theory. Therefore, even in this simplified example, one cannot establish any link between 
the sign of the Universe's baryon-number asymmetry with the observable CP-violating 
phases at low energies. 

In the very special case where hij has two zeros, one has only one CP-violating 
phase. In this case the CP-violating phase in neutrino oscillations is connected with the 
phase in Leptogenesis or, equivalently, the sign of the baryon-number asymmetry in the 
Universe [89]. 16 Thus, in this case, one may indeed test directly the idea of Leptogenesis. 
It is interesting that such a restricted model, where h\% — ^21 = 0, is still consistent with 
data on neutrino oscillations. 

5 Dark Matter Considerations 

It is certainly possible that the mechanism that generates a primordial matter-antimatter 
asymmetry in the Universe is not physically related to the existence of a non-luminous 
component of the energy density of the Universe, a component that now accounts for 
about 25 % of the total energy density. However, it would be very interesting if these 
two phenomena, so central to the history of the Universe, were connected in some deep 
way. It turns out, as we shall see, that if Leptogenesis is the mechanism by which a 
primordial matter-antimatter asymmetry in the Universe is established, it considerably 
impacts what the dark matter in the Universe can be. 

Of the three viable options for dark matter, from the point of view of particle physics, 

16 The prediction of the sign of the CP-violating phase in neutrino oscillations depends on which heavy 
Majorana neutrino is responsible for Leptogenesis. This problem is solved in the inflaton-decay scenario 
in supersymmetry (SUSY) theories, because one choice is unable to produce enough lepton-number 
asymmetry due to the constraint on the reheating temperature Tr < 10 7 GeV [90]. 
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two are either linked or constrained by thermal Leptogenesis and the third has clear 
connections to nonthermal Leptogenesis. Before discussing these points in some detail, 
it is useful to briefly review the extant dark matter candidates motivated by particle 
physics. 

5.1 PQ Symmetry and Axions 

It is well known that QCD admits the presence of an additional CP violating term in its 
Lagrangian density [17], 



where F iiiv = l/2e a p^ u F^ . If 9 is non-vanishing C g , which is C-even and P-odd, vi- 
olates CP and T invariance. Because possible CP violating parameters of the strong 
interactions, like the electric dipole moment of the neutron, are very tightly bounded by 
experiment [91], the parameter 9 must be very small (9 < 1CT 10 ) [92]. The reason for 
this is a mystery, and is known as the Strong CP Problem [93]. 

Probably the most 'natural' solution suggested for the Strong CP Problem is to 
assume that the total Lagrangian for the strong and electroweak interactions is invariant 
under a global chiral U (1) pq symmetry [94]. Even though this symmetry is spontaneously 
broken, one can show [94] that as a result of the U(1)pq symmetry the parameter 9 is 
driven to zero. In effect, what happens is that the CP violating Lagrangian term Cq is 
replaced by a CP conserving interaction between the CP odd pseudo-Goldstone boson 17 
associated with the spontaneous breakdown of U(1)pq - the axion - [95] and FF: 



Originally, it was supposed [94] that the U(1)pq symmetry was broken at the elec- 
troweak scale. Then f a ~ vf — 200 GeV and the axion mass lies in the keV range. 
However, axions in this mass range, which are coupled with strength l/f a , have been 
ruled out by experiment [93] . Astrophysical considerations, however, impose very strong 
constraints on axions much lighter than a keV, as their emission from stars would sig- 
nificantly alter their properties. Only axions that are sufficiently weakly coupled (hence, 
with large enough f a and thus a correspondingly small mass) avoid these constraints, 
and one finds the bound [96] f a > 10 10 GeV. On the other hand, f a cannot be arbitrar- 
ily large, because zero-momentum axion oscillations in the early Universe would carry 

17 Axions are not true Goldstone bosons because the U(1)pq symmetry is anomalous [95]. In fact, the 
same effective potential for axions that serves to drive 9 to zero gives axions a small mass. This mass 
is of order [93] rr? a ~ m q h^ CD / f%, where f a is the scale where the U(1)pq symmetry breaks down 
spontaneously, and m q is the (light) quark mass. 




(90) 



8tt/. ' ^ 
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enough energy density (proportional, approximately, to f a ) to overdose the Universe [97]. 
Thus, for an appropriate value for f a , axions can be the dark matter in the Universe. In 
particular, one finds [98] that f a ~ 10 12 GeV gives tt a ~ 1. 

5.2 Dark Matter Candidates from Super symmetry 

Supersymmetry, a boson-fermion symmetry, has been invoked extensively as the solution 
of the so-called hierarchy problem. This problem is related to the fact that without some 
stabilizing influence radiative corrections in the electroweak theory would naturally push 
the Fermi scale vf to have the value of whatever cutoff delimits the validity of the theory. 
Typically, this cutoff is imagined to be at the Planck scale Mp, and why vp « Mp is the 
hierarchy problem. This problem is resolved if there is some low energy (spontaneously 
broken) supersymmetry. Due to the fermion-boson nature of supersymmetry, radiative 
corrections of parameters in the electroweak theory (like vp) are now only logarithmically 
dependent on the cutoff, not quadratically dependent. Hence, effectively, if there is some 
low energy supersymmetry one can contemplate having a hierarchy like Vf « Mp, 
because radiative shifts can only change vp logarithmically. 

In general, supersymmetric theories possess a discrete symmetry (R-symmetry) that 
distinguishes particles from their supersymmetric partners. As a result, the lightest su- 
persymmetric particle (the LSP) is stable and, in principle, could be the source of the 
dark matter in the Universe. Indeed, it is known [64] that the energy density of particles 
of mass of 0(vf), whose annihilation cross section is of electroweak strength, is of the 
order of the critical energy density that closes the Universe. With supersymmetric part- 
ners of ordinary particles having electroweak scale masses and interactions, the LSP is 
therefore an ideal candidate for the dark matter in the Universe [99]. In this review we 
will discuss both the cases of neutralinos (the SUSY partners of gauge and Higgs boson) 
and of gravitinos (the spin 3/2 partner of the graviton) as LSP candidates. 

5.3 Extended Structures 

Scalar fields are necessary ingredients of the standard electroweak model, as well as its 
supersymmetric extension. It is well known that theories with scalar fields can lead to the 
formation of nontopological solitons. These extended structures, known as Q-balls [100], 
may be stable or unstable and arise when some scalar field carries a conserved U{\) 
charge. For example, in supersymmetric theories sleptons and squarks carry, respectively, 
lepton and baryon number. 

In supersymmetric theories, more generally, Q-balls can develop along flat directions 
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of the scalar potential [101]. These Q-balls can, in a number of instances, carry baryon 
number. If the baryon number of the Q-ball is large enough, and its mass is small enough, 
the baryonic Q-balls are stable. Because of their stability, one can imagine that these 
Q-balls could be the dark matter in the Universe. 18 Typically [101], if stable Q-balls exist 
they have both very large baryon number (B ~ 10 26 ) and are very massive (Mq < 10 26 
GeV). Unfortunately, this makes their detection very difficult, because their flux is very 
low [104]. 



5.4 Natural Connection of Axions with Leptogenesis 

The scale of U(1)pq breaking needed for axions to be the dark matter in the Universe 
(/„ ~ 0.3 x 10 12 GeV) is close enough to the mass of the lightest right handed neutrino 
(Mi ~ 10 10 GeV) needed for Leptogenesis to seek for a common linkage. In fact, the 
existence of such a linkage was observed long ago by Langacker, Peccei and Yanagida 
[105]. What these authors observed was that if M\ were due to the VEV of a scalar field 
a, one could identify this field as carrying a PQ-symmetry rather that lepton number. 

Let us examine this assertion in a bit more detail by looking at the Yukawa interac- 
tions of the quarks and leptons with the three Higgs fields 19 fa, fa and a. Schematically, 
one has 

^Yukawa = h a aN R N R + hN R fa£ L + h u u R faq L + fdd R faq L + fe R fa£ L + h.c. . (92) 

One sees that Eq. (|9"2*)) is invariant under a PQ-symmetry, where 

fa, fa -> e ia fa,e ia fa (93) 

N R ,l R ,u R ,d R -> e ia N R ,e ia l R ,e ia u R ,e ia d R , (94) 

provided that 

a -> e- 2ia a. (95) 

To allow < a >= f a >> Vp, as in all invisible axion models [106,107], requires one 
fine tuning. In the above case, this requires the PQ-invariant term in the scalar potential 

V = Kofafa + h.c. (96) 

18 This, however, is not easily achieved because, in general, the squarks are unstable with their baryon 
number eventually residing on quarks. If the squarks are light enough, stability can be achieved. However, 
as Kasuya et al. point out [102], it is difficult to explain both the baryon asymmetry and the dark matter 
density simultaneously. Nevertheless, there are scenarios where unstable Q-balls are the source for both 
baryogenesis and neutralino dark matter [103]. 

19 For a PQ symmetry to exist one needs to have two SU(2) doublet Higgs fields, 4>\ and </>2, rather 
than just the single Higgs field of the Standard Model H (and its Hermitian adjoint W). 
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to have the constant k ~ v F /f a , to allow electroweak symmetry breaking to occur at a 
scale much below the scale of U(1)pq symmetry breaking [vp « f a )- 

The Yukawa couplings of this model guarantee that the mass of the lightest right- 
handed neutrino and f a are related: Mi = 2(h a ) u f a . Thus, if axions are the source of 
the dark matter energy density in the Universe and the baryon aymmetry arises from 
Leptogenesis, because Odm ~ fa and Qb ~ Mi ~ f a , their ratio is independent of the 
scale of Upq(1) breaking. Hence it is perhaps not surprising that this ratio is of order 
unity. 



5.5 The Gravitino Problem in Supersymmetric Theories 

As we discussed in Section 4, for Leptogenesis to be effective, the mass of the lightest 
right-handed neutrino has to be greater than 2 x 10 9 GeV. This bound, in turn, means 
that thermal Leptogenesis must have occurred at temperatures above 2x 10 9 GeV. Hence, 
if the Universe went through an inflationary period, as all evidence seems to suggest [4] , 
the reheating temperature after inflation Tp must have been greater than 2 x 10 9 GeV for 
Leptogenesis to be the source of the matter-antimatter asymmetry in the Universe. This 
high reheating temperature is problematic for supersymmetric theories because it leads 
to an overproduction of light states, like the gravitino, with catastrophic consequences for 
the evolution of the Universe after inflation. Unless these observational inconsistencies 
can be avoided, it appears that Leptogenesis in supersymmetric theories cannot produce 
the desired baryon asymmetry in the Universe. 

The production of gravitinos after inflation has been studied in some detail [108] . The 
thermal production of gravitinos produced by the strong interactions of quarks, squarks, 
gluons and gluinos is governed by the Boltzmann equation [109] 

dn 3/2 



dt 



+ 3Hn 3/2 = C 3/2 (T), (97) 



where 

r m _ 3C(3)a 3 (T) T 6 ( roj \ 

° 3/2{T) - — ^ — w P y 1 + 3^ J F{T) - (98) 

Here F(T) is a thermal factor of O(10) and rrig and m 3 / 2 are, respectively, the gluino and 
the gravitino masses. Integrating Eq. (|9*7j) to a reheating temperature T R , the resulting 
relic density of produced gravitinos is given by 



n 3/2 h 2 ~ 0A4a 3 (T R ) 




Tr \ ( m ^ \ (99) 
10 10 GeW UOOGeW ' 1 ; 



where h is the scaled Hubble parameter and /i ~ M z . 
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If gravitinos are stable (i.e. they are the LSP), the WMAP constraint on the amount 
of dark matter in the Universe [4] 

n DM h 2 = 0.11261°;™* (100) 

constrains VLz^h 2 to be below this value and, for any given reheating temperature Tr 
and gravitino mass m 3 / 2 , gives a bound on the gluino mass. If, on the other hand, the 
gravitinos are not stable, their rate of production for Tr > 2 x 10 9 GeV is so large 
that subsequent gravitino decays completely alter the standard Big Bang Nucleosyn- 
thesis (BBN) scenario. Thus, in either case, there are severe constraints imposed on 
supersymmetric dark matter, which we will discuss in detail below. 

If the gravitino is unstable, it has a long lifetime and decays during or after BBN 
for an interesting range of the gravitino mass, 777,3/2 — 100 GeV — 10 TeV. The gravitino 
decay products destroy the light elements produced by the BBN and hence the relic 
abundance of gravitinos is constrained from above to keep the success of the BBN [110]. 
This leads to an upper bound of the reheating temperature Tr after inflation, since the 
abundance of gravitinos is proportional to the reheating temperature. A recent detailed 
analysis derived a strigent upper bound Tr < 10 6 ~ 7 GeV when the gravitino decay has 
hadronic modes (see Fig. (J3J)) [90]. This upper bound is much lower than the temperature 
for Leptogenesis, T R > 2 x 10 9 GeV [69, 72]. Therefore, thermal Leptogenesis seems 
difficult to reconcile with low energy supersymmetry if gravitino masses lie in the range 
7773/2 — 100 GeV — 10 TeV - a natural range for Supergravity (SUGRA) models. 

5.6 Solutions to the Gravitino Problem in Thermal 
Leptogenesis 

There have been several attempts to solve the gravitino problem in thermal Leptogenesis. 
Here we will briefly review a number of these proposed solutions. 

One possibility has been proposed by Pilaftsis who considers quasi-degenerate heavy 
Majorana neutrinos (Mi ~ M2) [111]. In this model the lepton-asymmetry parameter e 
is enhanced by a factor of Mi /(Mi — M 2 ) and hence the decays of both Ni and N 2 may 
produce enough asymmetry even for Tr < 10 6 ~ 7 GeV. However, it is difficult to find a 
compelling justification for having such a degeneracy in the heavy neutrino spectrum. 

Another proposal was made by Bolz, Buchmuller and Plumacher [112], who consider 
the case where the gravitino is the stable lightest SUSY particle (LSP). In this case 
the next lightest supersymmetric particle (NLSP) is the subject of the cosmological 
constraint, because its decay products may destroy the light elements created by the 
BBN, much like the unstable gravitino. Detailed analyses show that this scenario favors 

36 



lo glo( T 3/ 2 / SeC ) 

4 2 




10 J 



iL 



95% C.L. 

77 = (6. 1 ±0.3)x 10" 10 _ 



I I 



io £ 



io 3 

IT1 



3/2 



10 4 
(GeV) 



io f 



IO 10 

IO 9 

IO 8 

10 v 

IO 6 

10 5 

IO 4 

10 3 

IO 2 



Figure 5: Upper bounds on the reheating temperature as function of the gravitino mass 
for the case where the gravitino dominantly decays into a gluon-gluino pair. From [90]. 



the f NLSP compared to the neutralino NLSP and gravitino masses below m 3 / 2 — 
100 GeV [113, 114]. 20 In general, the gravitino production can be dominated by NLSP 
decays [115] or by thermal processes [116]. 

A third proposed solution makes use of gauge-mediation model of supersymmetry 
breaking in which the gravitino is the stable LSP with a mass 7713/2 < 1 GeV. It turns 
out that, if the gravitino mass is 7773/2 < 16 eV [117], then there is no gravitino problem. 
However, in this case the gravitino cannot be the dark matter. It must be something else, 
perhaps the axion. For the range of gravitino masses T77 3 / 2 ~ 100 keV — 1 GeV, there is 
the interesting possibility that late-time entropy production in a class of gauge mediation 
models can naturally make the gravitino the dominant component of dark matter [118]. 
In this scenario the reheating temperature can be as high as Tr ~ 10 13 GeV. A light 
axino with mass (9(1 keV) as LSP and a gravitino as NLSP would solve the gravitino 
problem [119]. 

Finally, another possible solution arises if supersymmetry breaking effects are trans- 
mitted to the Standard Model sector through the scale anomaly, resulting in very heavy 
gravitino masses m 3 / 2 > 100 TeV. In this case the gravitino decays before the time of 
BBN and hence there is no cosmological problem. However, the gravitino decay modes 
contain always one LSP and hence the relic abundance of the gravitino must be con- 

20 The gravitino mass is even less cnstrained if the LSP is a scalar neutrino or the gluino. For a 
class of supergravity models an upper bound of 5 x 10 9 GeV on the reheating temperature has been 
obtained [114]. 
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strained from above so that the density of the nonthermal LSP produced by the grav- 
itino decays does not exceed the dark-matter density. This condition leads to T R < 10 11 
GeV [120,121], which is consistent with thermal Leptogenesis. 

We stress here that each of the above 'solutions' predicts distinct particle spectra 
at the TeV scale, which may be testable in future collider experiments at the Large 
Hadron Collider (LHC). If one discovers supersymmetry but all the above possibilities 
are excluded experimentally, this would argue strongly against thermal Leptogenesis 
although nonthermal Leptogenesis could still be viable. On the other hand, if some of 
these scenarios are confirmed experimentally, thermal Leptogenesis will become much 
more compelling. 

5.7 Leptogenesis and Lepton Flavor Violation in SUSY Models 

Lepton flavor violation is another area in which thermal Leptogenesis and supersymme- 
try may have some linkages. If the neutrinos have a mass, lepton flavor violation (LFV) 
processes such as /i — > e + 7 decay can occur. In non supersymmetric theories, these 
processes are strongly suppressed by a factor of (m u / Mw) 2 in rate and hence are un- 
measurable physically. However, this is not the case in the SUSY Standard Model, if the 
seesaw mechanism is effective. 

In the SUSY Standard Model scalar quarks and leptons are assumed to have a uni- 
versal SUSY-breaking soft mass, m , at the Planck scale. Otherwise one would have too 
large flavor-changing neutral currents (FCNC). However, even then quantum corrections 
resulting from Yukawa interactions of the quarks generate a violation of the universal- 
ity of the soft masses for scalar quarks, which induces FCNC. In the lepton sector the 
Yukawa couplings of the superheavy Majorana neutrinos iVj also generates non-universal 
masses for scalar leptons that serves as a source of LFV [122]. 21 If the relevant Yukawa 
couplings are of 0(1), or equivalently M3 ~ 10 15 GeV, one predicts a branching ratio 
for n — > e + 7 decay that may be testable in future experiments. However, an accurate 
prediction for LFV processes is very difficult, since it hinges on unkown Yukawa coupling 
constants [124]. In particular, the constraint coming from Leptogenesis that M\ ~ 10 10 
GeV is not strong enough to suggest that LFV processes have potentially testable rates. 

21 The Yukawa couplings hij of the Higgs to the A^s and leptons induce flavor dependent 
soft masses for scalar leptons. At the one-loop level the induced mass is given by (m^)fj ~ 
— 6m§/ \A^) 2 h\ k hkjln{Mp/Mk), where too is the universal soft mass for scalar leptons at the Planck 
scale Mp. Thus, one may obtain information on the high-energy Yukawa coupling hij and the heavy 
neutrino masses Mk by measuring directly the mass matrix for the scalar leptons. However, the phases 
in fvh are different from the phases contributing to Leptogenesis [123]. 
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6 Nonthermal Leptogenesis 



Supersymmetry is an important symmetry for the unification of all interactions and all 
matter, and the SUSY Standard Model is considered as a plausible scenario for producing 
new physics at the TeV scale. Thus, it is quite interesting to consider theories where 
supersymmetry is spontaneously broken in a hidden sector connected to ordinary matter 
by gravitational strength interactions- the SUGRA framework. The seesaw mechanism 
is easily incorporated into this framework. However, as we discussed in some detail in 
Section 5, the gravitino problem argues against thermal Leptogenesis, particularly in 
SUGRA. 

A possible solution to this problem may be provided by nonthermal Leptogenesis 
[10,125-129], where one does not have a strong constraint on the reheating temperature. 
We will discuss here specifically nonthermal Leptogenesis via inflaton decay [125, 126], 
which we consider an interesting scenario. In the next subsection, we present general 
arguments for this scenario and show that it suggests a lower bound on the mass of the 
heaviest light neutrino m 3 > 0.01 eV. In the subsequent subsection, we will also discuss 
the Affleck-Dine mechanism [10] for Leptogenesis which, specifically in supersymmetric 
theories, is also an interesting mechanism to generate the matter-antimatter asymmetry. 



6.1 Nonthermal Leptogenesis via Inflaton Decay 

Inflation early on in the history of the Universe is one of the most attractive hypothesis 
in modern cosmology, because it not only solves long-standing problems in cosmology, 
like the horizon and the flatness problems [130], but also accounts for the origin of 
density fluctuations [131]. In this subsection we discuss the hypothesis that the inflaton 
$ decays dominantly into a pair of the lightest heavy Majorana neutrinos, $ — ► Aq + Aq. 
We assume, for simplicity, that other decay modes including those into pairs of N2 and A^ 
are energetically forbidden. The produced Aq neutrinos decay subsequently into H + £l 
or W + £ L . If the reheating temperature Tr is lower than the mass Mi of the heavy 
neutrino Aq, then the out-of-equilibrium condition [11] is automatically satisfied. 

The above two channels for Aq decay have different branching ratios when CP con- 
servation is violated. Interference between tree-level and one-loop diagrams generates a 
lepton-number asymmetry [57]. Following our discussion in Section 4, the lepton asym- 
metry parameter e can be written as [128, 129, 132] 22 

3 M i t \ 

£ _= -^w m3eS ' (101) 

22 Because of supersymmetry, the asymmetry parameter e below is a factor of 2 larger than that given 
in Eq. JSHJ- 

39 



where the effective CP-violating phase S e s is given by 
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(102) 



Numerically, one obtains for the e parameter 



The chain decays $ — ► i\q + iVi and N\ —>■ H + or iff + ^ reheat the Universe 
producing not only the lepton-number asymmetry but also entropy for the thermal bath. 
The ratio of the lepton number to entropy density after reheating is estimated to be [126] 

ul_ _ 3^ T R 
s 2 m$ 

10 (lo^)©(aSv)> <-> 

where m$ is the inflaton mass and we have taken 5 e g = 1. This lepton-number asymmetry 
is converted into a baryon-number asymmetry through the sphaleron effects and one 
obtains [35] 

T " "2^T ■ (105) 

We should stress, here, an important merit of the inflaton-decay scenario: It does 
not require a reheating temperature Tr ~ Mi, but it requires only m$ > 2Mi. On 
the other hand, for thermal Leptogenesis to work it is necessary that Tr ~ Mi, which 
necessitates higher reheating temperature for Leptogenesis to produce enough matter- 
antimatter asymmetry. 

If one assumes that Tr < 10 7 GeV to satisfy the cosmological constraint on the 
gravitino abundance [90] discussed earlier and uses m$ > 2Mi, the observed baryon 
number to entropy ratio [4] gives a constraint on the heaviest light neutrino: 

m 3 > 0.01 eV. (106) 

It is very interesting that the neutrino mass suggested by atmospheric neutrino oscillation 



experiments, JAml tra ~ 0.05 eV, just satisfies the above constraint. However, to get 
this bound we assumed that the inflaton decays dominantly into a pair of iV\s. If this 
branching ratio is only 10 %, the lower bound on the neutrino mass exceeds the observed 



neutrino mass y Ami tm ~ 0.05 eV. 

A variety of models have been considered to restore the bound of Eq. (|106J) by 
imposing a symmetry. However, it is perhaps most interesting to consider that the scalar 
partner of the heavy Majorana neutrino N\ is the inflaton itself [127], and the inflaton 
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decay into a lepton plus a Higgs boson gives an effective branching ration of 100%. In 
this model, one must assume that the initial value of the scalar partner of Ni is much 
larger than the Planck scale to cause inflation (chaotic inflation [133]). However, chaotic 
inflation is not easily realized in SUGRA, because the minimal supergravity potential 
has an exponential factor, exp(0*0/M^), that prevents any scalar field <fi from having 
a value larger than the reduced Planck scale Mq — 2.4 x 10 18 GeV. Ref. [134] uses a 
restricted form of the Kahler potential. 

6.2 Affleck-Dine Leptogenesis 

In the SUSY Standard Model, in the limit of unbroken supersymmetry, some combina- 
tions of scalar fields do not enter the potential, constituting so-called flat directions of 
the potential. Since the potential is almost independent of these fields, they may have 
large initial values in the early Universe. Such flat directions receive soft masses in the 
SUSY-breaking vacuum. When the expansion rate H cxp of the Universe becomes com- 
parable to their masses, the flat directions begin to oscillate around the minimum of the 
potential. If the flat directions are made of scalar quarks and carry baryon number, the 
baryon-number asymmetry can be created during these coherent oscillations. This is the 
Affleck-Dine (AD) mechanism for Baryogenesis [10]. 

QCD corrections, however, make the potential of the AD fields milder than |0| 2 . This 
allows non-topological soliton solutions (Q-balls) [135] to form in the early Universe, as 
a result of the coherent oscillations in the flat directions. Because Q-balls have long 
lifetimes, their decays produce a huge amount of entropy at late times. To avoid this 
problem one must choose parameters in the SUSY theory so that the density of the 
lightest SUSY particle (LSP) does not exceed the dark matter density in the present 
Universe [135]. Although this may not be a problem, it is much safer to consider flat 
directions without QCD interactions, because such directions most likely do not have 
Q-ball solutions. 

The most interesting candidate [128] for such a flat direction is 

0, = (2^) 1 / 2 , (107) 

where £i is the lepton doublet field of the i-th family. Here, H and l x represent the scalar 
components of the corresponding chiral multiplets. The Yukawa interactions of H make 
the potential of (pi steeper than the mass term and hence there is no instability of the 
coherent oscillation (i.e. there are no Q-ball solutions). Because this flat direction carries 
lepton number, a lepton asymmetry will be created during the coherent oscillation (AD 
Leptogenesis) [128]. Sphaleron processes then transmute, in the usual fashion, this lepton 
asymmetry into a baryon asymmetry. 
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The seesaw mechanism induces a dimension-five operator in the superpotential for 
the theory, 23 

W = ^rMHf, (108) 



where we have used a basis in which the neutrino mass matrix is diagonal. With this 
superpotential we have a SUSY-invariant potential for the flat direction given by 

2 

Tfl 

^susy = ^y^l 6 - (109) 
In addition to the SUSY-invariant potential we have a SUSY-breaking potential, 

SV = ml\<Pn^^(a m ^ + h. C .). (110) 

Here, a m is a complex number. We take ~ msusY — 1 TeV and \a m \ ~ 1. The second 
term in SV is very important, because it gives rise to the lepton- number generation. 

We assume that the flat direction acquires a negative (mass) 2 induced by the 
inflaton potential and rolls down to the point balanced by the SUSY-invariant potential 



Vsusy during inflation. Thus, the AD field has an initial value of y-ffmf | (H)\ 2 /m u , where 
ifinf is the Hubble constant (the expansion rate) during inflation. decreases in amplitude 
gradually after inflation, and begins to oscillate around the potential minimum when the 
Hubble constant H exp of the Universe becomes comparable to the SUSY-breaking mass 
m^. At the beginning of the oscillation, the AD field has the value \<f> \ ~ ^Jm^\(H)\ 2 /m u 
which, as shown below, is an effective initial value for Leptogenesis. 

Let us consider now lepton-number generation in this scenario. The evolution of the 
AD field is described by 

^ + 3iJ cxp - + — -0, (111) 
where V = Vsusy + SV. Because the lepton number is given by 

n L = ,(^0-0^j , (112) 

the evolution of is given by 

dn L m S u SY m u . 

— + 3H cxp n L = — — ^lm(a m ) . (113) 

The motion of in the phase direction generates the lepton number. This is pre- 
dominantly created just after the AD field starts its coherent oscillation, at a time 



23 



For ease of notation we have dropped the subscript i below. 
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t osc ~ l/H osc ~ because the amplitude \(f>\ damps as t 1 during the oscillation. 

Thus, we obtain for the lepton number 

n L - ol/u\l2 d eS\a m <Po\ X t osc , (114) 

2 \{ H )\ 

where 8 e g = sin(4arg0 + arga m ) represents an effective CP-violating phase. Using 



m S usY - rut, |0 O | - \J m 4>\{ H )\ 2 / m v and t osc ~ 1/m^, we find 

n L * <W^J^ . (115) 

After the end of inflation, the inflaton begins to oscillate around the potential mini- 
mum and nL/Pinf stays constant until the inflaton decays. Here p; n f is the energy density 
of the inflaton. The inflaton decay reheats the Universe producing entropy s. Because 
p/s = 2Tr/4, we find for the lepton-number asymmetry the expression 

s ~\s ) \ Pmf J ~ 4:M G 6m u M G ^ LLb) 

Here we have used p in{ ~ 3m^M^- at the begining of the AD field oscillation (when 
most of the lepton number is generated). This lepton-number asymmetry is converted 
to a baryon-number asymmetry by the KRS mechanism. In this way one obtains for the 
baryon-number asymmetry 

n B ^ 1 \{H)\ 2 T R 



s 23 m v M 2 G 



(117) 



The observed ratio n^/s ~ 0.9 x 1CT 10 implies m u ~ 10~ 9 eV for Tr ~ 10 6 GeV. This 
small mass corresponds to the mass of the lightest neutrino. We should note that for 
such a low reheating temperature one may neglect the effects due to thermal mass for 
the AD field [136]. 



7 Conclusions and Summary of Results 

In this article we have discussed the physical mechanism responsible for the origin of 
matter in the Universe. Both the rather large observed value for the ratio of baryons to 
photons, t]b, in the present epoch and the absence of antimatter are the consequences of a 
primordial asymmetry between matter and antimatter generated early on in the Universe. 
Although a variety of mechanisms have been proposed for producing this primordial 
asymmetry, in this review we have focused on Leptogenesis as the origin of matter. In 
our view, this is the most appealing scenario for the origin of matter, for at least three 
reasons: 
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1) Explicit lepton number violation is very natural once one includes right-handed 
neutrinos in the Standard Model. Furthermore, the lightness of the observed neutrinos 
strongly suggests, through the seesaw mechanism, the presence of superheavy neutrinos, 
whose decays can produce a lepton-antilepton asymmetry. 

2) Quantum mechanically, through the KRS mechanism, one can automatically turn 
a leptonic asymmetry into a baryonic asymmetry. Indeed, because of the existence of 
these sphaleron processes, the origin of matter is linked to phenomena in the early Uni- 
verse that result in the establishment of a (B-L)-asymmetry, like Leptogenesis. 

3) If neutrino masses lie in the range 10~ 3 eV < im < O.leV, as suggested by neutrino 
oscillation experiments, the leptonic asymmetry produced in thermal Leptogenesis is both 
independent of the abundance of heavy neutrinos and of any pre-existing asymmetry and 
has the right magnitude to yield the observed value for r\ B . 

Because, in the final analysis, t)b is just one number, it is important to ask if the 
particular mechanism proposed for the origin of matter has other consequences. Thus 
in this review we examined in some detail how Leptogenesis fit with ideas proposed to 
explain the dark matter that constitutes about 25% of the Universe's energy density. 24 
We pointed out that axionic dark matter is perfectly compatible with Leptogenesis. 
Indeed, it is possible to very naturally link the scale of the heavy neutrinos with that of 
U{\)pq breaking f a , so that the ratio f2 DM /f2 B is independent of these large scales. The 
situation, however, is more complex in the case of supersymmetric dark matter. 

To be effective, thermal Leptogenesis needs to occur at high temperatures, above 
T = 2 x 10 9 GeV. This means that the Universe after inflation must have reheated to 
at least this temperature. However, in supersymmetric theories such a high reheating 
temperature is problematic as it leads to an overproduction of gravitinos. When they 
decay, gravitinos of such abundances completely alter the primordial abundance of light 
elements produced in Big Bang Nucleosynthesis. The gravitino problem, however, is 
not fatal as there are a number of ways to mitigate the overproduction of gravitinos. 
Nevertheless, if Leptogenesis is at the root of the origin of matter, the supersymmetric 
spectrum at low energies and the nature of the LSP are quite constrained. Thus, in a 
sense, Leptogenesis is also quite predictive in this context. 

Although much of our review, very naturally, focused on thermal Leptogenesis, we 
also discussed two examples where matter originated through a leptonic asymmetry pro- 
duced in nonthermal processes. These models, although much more speculative, illustrate 
some of the possible other options for the origin of matter. Naturally, in this case some 
of the specific predictivity is lost. 

24 Wc did not try to examine models of dark energy in the light of Leptogenesis, because our under- 
standing of dark energy is still in its infancy. 
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